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Abstract. We analyse the 267 radio sources from our deep (flux limit of 
42 /iJy at the field center at 1.4 GHz) Chandra Deep Field South 1.4 and 5 
GHz VLA survey. The radio population is studied by using a wealth of multi- 
wavelength information, including morphology and spectral types, in the radio, 
optical, and X-ray bands. The availability of redshifts for ~ 70% of our sources 
allows us to derive reliable luminosity estimates for the majority of the objects. 
Contrary to some previous results, we find that star-forming galaxies make up 
only a minority (ss 1/3) of sub-mjy sources, the bulk of which are faint radio 
galaxies, mostly of the Fanaroff-Riley I type. 



1. The VLA-CDFS Survey 



The extragalactic radio source population ranges from normal galaxies with lu- 
minosities ~ 10 20 W/Hz to galaxies whose radio emission is as much as 10 4-7 
times greater owing to regions of massive star formation or to an active galactic 
nucleus (AGN). The population of radio sources in the sky with flux densities 
> 1 mJy is dominated by AGN driven emission generated from the gravita- 
tional potential associated with a supermassive black-hole in the nucleus. For 
these sources, the observed radio emission includes the classical extended jet 
and double lobe radio source as well as compact radio components more directly 
associated with the energy generation and collimation near the central engine of 
the AGN. Below 1 mJy there is an increasing contribution to the radio source 
population from synchrotron emission resulting from relativistic plasma ejected 
from supernovae associated with massive star formation in g alaxies or groups 
of ga l axies, often associated with mergers or inte ractions (e.g.. IWindhorst et al.l 
19951 : iRichards et alj[l998l ; iFomalont et al.|[200l ). However, the mix of star for- 
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mation and AGN related radio emission and the dependence on epoch is not 
well determined, although the commonly accepted view ascribes to star- forming 
galaxies a major role at these faint fluxes. 

The Chandra Deep Field South (CDFS) area is part of the Great Observa- 
tories Origins Deep Survey (GOODS) and as such is one of the most intensely 
studied re gion of the sky. High sensitivity X-ray ob servations are available from 
Chandra (iGiacconi et al.ll2002l : ITehmer et al.ll2005h . The GOODS and Galaxy 
Evolution from Morphology and Spectral energy distributions (GEMS) multi- 
band imaging programs using the HST A dvanced Camera for Surveys (ACS) giy e 
sensitive high resolution optical images (|Giavalisco et al.ll2004l : [Rix et alJl2004l ). 
Ground based imaging and spectroscopy are available from the ESO 2.2m and 
8m telescopes, and infrared observations from the Spitzer space Telescope. 

We observed the CDFS with the NRAO Very Large Array (VLA) for 50 
hours at 1.4 GHz mostly in the BnA configurations in 1999 and 2001, and for 32 
hours at 5 GHz mostly in the C and CnB configurations in 2001. The effective 
angular resolution was 3.5" and the minimum root-mean-square (RMS) noise as 
low as 8.5 //Jy per beam at 1.4 GHz and 7 /xJy per beam at 5 GHz. A total of 
267 radio sources were catalogued at 1.4 GHz, 198 of which are in a complete 
sample with signal-to-noise ratio (SNR) greater than 5 and located within 15' 
of the field center. The corresponding flux density limit ranges from 42 //Jy at 
the field center to 125 fiJy near the field edge. These deep radio observations 
complement the larger area, but less sensit ive lower resolution observations of 
the CDFS discussed bv lAfonso etaD (|2006l ). 

Our set of ancillary data is quite unique and allows us to shed new light 
on the nature of the sub-mJy radio source population. These data include: a) 
reliable optical/near-IR identifications for ~ 96% of the radio sources; b) optical 
morphological classification for ~ 54% of the sample; of these: 38% are spiral 
or interacting, 30% are elliptical or lenticular, 16% are compact, and 15% are 
irregular or complex; c) redshift information for 186 (~ 70%) of the objects: 
85 spectroscopic (32%) and 101 photometric (38%). The redshift distribution 
covers the range 0.03 - 3.7 and peaks at z ~ 0.6 — 0.7, with a sharp cut-off at 
z ^> 1. The mean redshift is {z) ~ 0.8; d) X-ray detections for 85 (~ 32%) of 
the objects, upper limits for all the others. 

We present here preliminary results on the sub-mJy source population of 
the VLA-CDFS. More details on the VLA-CDFS can be found in a series of four 
papers, which address the radio data (Kellermann et al. 2007, in preparation), 
the optical and near IR counterparts to the observed radio sources (Mainieri et 
al. 2007, in preparation), their X-ray spectral properties (Tozzi et al. 2007, in 
preparation), and the source population (Padovani et al. 2007, in preparation). 
Throughout this paper spectral indices are written S u oc v~ OL and the values 
Hq = 70 km s _1 Mpc -1 , Qm = 0.3, and Q\ = 0.7 have been used. 



2. The sub-mJy number counts 

Fig. Q] shows the counts of sources at 1.4 GHz over the entire range of observed 
flux densities. The highest density points are based on a compilation from 
Condon (2006, p. a). Various studies over many decades have shown that there 
are three regions of the counts: a) > 1 Jy, where the counts are rising with respect 
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to Euclidean counts, due to the strong evolution of luminous quasars and radio 
galaxies; b) 3 mJy - 1 Jy, still dominated by quasars and radio galaxies, and 
characterized by a drop-off predominantly caused by the redshift effects and the 
cutoff above z = 3 in the density of these luminous radio sources; c) < 3 mJy, 
where the counts flatten closer to Euclidean, well above the extrapolation from 
higher flux densities, and the field to field scatter appears to increase below 
about 400 /xJy This flattening of the slope of the counts, now observed in the 
flux density range 35 — 300^Jy, is produced by a population of objects, which 
is now dominant in the counts. The characterization of this population is the 
major goal in the study of weak radio sources and of our project. 
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Figure 1. The number counts at 1.4 GHz normalized to a Euclidean value. 
Our data points are shown as triangles. The dashed line, which is the fit to 
the CDFS counts, is a reasonable fit to the average counts from all surveys at 
the lowest flux densities. 



3. The sub-mJy population 

To shed light on the sub-mJy population, we start by evaluating a measure of the 
radio-to-optical flux ratio, R = log(S'i.4GHz/>S'v)- About 40% of our sources have 
k-corrected iiQ values above the star-forming regime, R <^ 1.4, with some of our 
sources reaching extremely high values of R, equivalent to radio-to-optical ratios 
p> 10 4 — 10 5 . It has been suggested that the radio-to-optical flux ra tio cannot be 
used to discriminate betw een AGN and star- forming galaxies (e.g.. lAfonso et aT 



20051 : iBarger et al.ll2006l ). We find this to be the case at low (R <; 1.4) values, 



where there is a large overlap between the two classes. However, above this value 
practically only AGN are present amongst the classified sources. Note that a 



values have been k-corrected using templates typical for that class or, in the case of sources 
without morphological classification, assuming a template which was the average between el- 
lipticals and spirals 
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sample with a brighter magnitude cutoff would get very different results. For 
example, for V 22, only ~ 10% of the sources would have R values beyond 
the star-forming regime. This is an important p oint, which might explain some 
previous results (see also iGruppioni et al.lll999T ). 
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Figure 2. Top: B = log(Si.4GHz/ Sy) vs. the 1.4 GHz radio power for 
our sources. The horizontal (long-dashed) line at R = 1.4 and the verti- 
cal (short-dashed) line at logP r = 24 indicate the maximum value for star- 
forming galaxies and the approximate dividing lines between radio-loud and 
radio-quiet AGN. The mean values (and ler) for the sources without red- 
shift information assuming z = 1 and z — 2 are also shown. Bottom: The 
2—10 keV X-ray power vs. the 1.4 GHz radio power for sources with red- 
shift and X-ray information. The horizontal (long-dashed) line at L x = 10 42 
erg/s indicates the dividing line between AGN and normal galaxies. The 
dot-dashed line is the L x — P r r elationship for nearby star-forming galaxies 
(Ran alli. Comastri &: S etti 2003), converted to our cosmology, while the dot- 
ted lines denote the locus of radio galaxies. Arrows represent X-ray upper 
limits. Solid lines are the best fits for the various classes taking upper limits 
into account. 

Fig. [21 top, which plots R versus the 1.4 GHz radio power for our sample, 
shows that the majority (92%) of spiral and interacting galaxies fall in the region 
typical of star-forming galaxies, i.e., P r <J 10 24 W/Hz; most of them (90%) have 
also R < 1.4. Ellipticals and lenticulars span a very large range of radio powers, 
consistent with that of low-luminosity (e.g., Fanaroff-Riley type I) radio-galaxies, 
while all compact sources are consistent with being radio-loud AGN. Therefore, 
despite the substantial overlap at low powers between ellipticals and spirals, the 
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various types have quite distinct radio powers. Namely, (logP r ) = 22.65 ±0.10, 
23.46 ±0.19, 24.03 ±0.20, and 23.59 ±0.11 for spirals and interacting, ellipticals 
and lenticulars, compact, and "other" sources (i.e., without morphological clas- 
sification) respectively. The latter have the largest R values and a distribution 
of radio powers which is significantly different (P > 99.99%) from that of spi- 
ral and interacting galaxies, but similar to that of ellipticals. Furthermore, the 
objects without redshifts have quite faint optical magnitudes and are therefore 
very likely characterized by higher than average redshifts. For estimated z = 1 
or z = 2, the bulk of these sources ends up in the region typical of radio galaxies. 

Under the assumptions that: a) star- forming galaxies are identified with spi- 
rals and interacting galaxies; b) star-forming galaxies alone occupy the bottom- 
left quadrant of Fig. [21 top, their fraction would be 62^% of the objects with 
redshift information. However, only ~ 40% of these sources have a definite spiral 
morphology, and their maximum fraction (excluding ellipticals and lenticulars) 
is bound to be ^ 3/4. Moreover, as most of the objects without redshift in- 
formation have relatively large R values and therefore are unlikely star-forming 
galaxies, the total fraction of such objects in the sample is < 45%. 

Fig. [21 bottom, shows L x (2 -10 ke V) vs. P r for the sources with red- 
shift. Using ASURV (jLavallev et al. 19921 ) to take into account the upper lim- 
its on L x , we find a correlation between the two powers for the whole sample 
(L x oc P r °- 92 ±o.08 ; p > 99.99%). Moreover: a) the large majority (80%) of spiral 
and interacting galaxies have X-ray detections consistent with powers typical of 
galaxies and not of AGN (<^ 10 42 erg/s); the best fit correlation using ASURV 
overlaps almost exactly with the X-ray-radio power relationship for nearby star- 
forming galaxies (dot-dashed line); b) the opposite is true for compact sources, 
which all have X-ray powers > 10 erg/s and ^ 30 times larger than the other 
classes, at a given radio power; this is consistent with the idea that these sources 
are mostly radio- loud AGN (see also Fig. [21 top); c) ellipticals and lenticulars 
span a very large range of X-ray powers (similarly to what observed in the radio 
band), consistent with that of radio-galaxies; again, there is substantial overlap 
at low X-ray powers between ellipticals and spirals, with at least 55% of the 
former having L x < 10 42 erg/s; d) the "other" sources cover a very large range 
in radio and X-ray powers, with the best fit correlation mostly into the radio- 
galaxies area. Finally, as was the case for the radio powers, the various types 
have quite distinct X-ray powers. Namely, (logL x ) = 40.65 ±0.15, 41.34 ±0.21, 
43.40 ± 0.18, and 41.25 ± 0.30 for spirals and interacting, ellipticals and lentic- 
ulars, compact, and "other" sources respectively. The X-ray power distribution 
for the sources without morphological classification is different (P > 97%) from 
that of spiral galaxies but similar to that of ellipticals. 

Similarly to Fig. [21 top, we can define the bottom-left quadrant as that of 
star-forming galaxies. If star-forming galaxies alone occupy this quadrant, their 
fraction would be 38l|j% of the objects with redshift and X-ray information. 
However, only about half of these sources have a definite spiral morphology, 
and their maximum fraction (excluding ellipticals and lenticulars) is bound to 
be <^ 3/4. On the other hand, some of the sources with upper limits above 
L x ~ 10 42 erg/s might have X-ray powers below that value. Therefore, the total 
fraction of star-forming galaxies in the sub-sample with redshift information is 
likely to be « 30%. 
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Figure 3. The percentages of various morphological types as a function of 
radio flux density for the complete sample. 



Fig. [3] shows the percentage of various morphological types as a function of 
radio flux density for the complete sample. The fraction of spirals and interacting 
galaxies increases from ~ 10% at 1 mJy to ~ 20% at the survey limit, while 
that of ellipticals, lenticulars, and compact sources (i.e., the AGN, shown in the 
figure as 100 — %) decreases from ~ 60% to ~ 30%. Hence, sub-mJy star- forming 
galaxies down to ~ 40/iJy can make up at most 70% of the total. However, given 
the evidence presented above, a stronger upper limit is provided by the fraction 
of sources having R < 1.4 excluding the sources with elliptical and lenticular 
morphology, that is <^ 45%. Therefore, star- forming galaxies make up only 
a minority, between ~ 20 and 45% 1/3), of sub-mJy sources, the bulk of 
which are faint radio galaxies, mostly at relatively low powers and hence of the 
Fanaroff-Riley I type. 

The VLA is operated by NRAO, a facility of the NSF under cooperative 
agreement with Associated Universities, Inc. 



References 

Afonso, J., et al. 2005, ApJ, 624, 135 
Afonso, J., et al. 2006, AJ, 131, 1216 

Barger, A. J., Cowie, L. L., & Wang, W.-H., 2006, ApJ, in press ( |astro-ph/ 0609374) 
Fomalont, E. B., et al. 2002, AJ, 123, 2402 
Giacconi, R., et al. 2002, ApJS, 139, 369 
Giavalisco, M., et al. 2004, ApJ, 600, L93 

Gruppioni, C, Mignoli, M., & Zamorani, G. 1999, MNRAS, 304, 199 

Lavalley, M., Isobe, T., & Feigelson, E. 1992, ASP Conf. Ser. 25: Astronomical Data 

Analysis Software and Systems I, 25, 245 
Lehmer, B. D., et al. 2005, ApJS, 161, 21 
Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 
Richards, E. A., et al. 1998, AJ, 116, 1039 
Rix, H.-W., et al. 2004, ApJS, 152, 163 
Windhorst, R. A., et al. 1995, Nature, 375, 471 



